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Abstract 

We investigate the vacuum expectation value of the fermionic current induced by a mag- 
netic flux in a (2-|-l)-dimensional conical spacetime in the presence of a circular boundary. On 
the boundary the fermionic field obeys MIT bag boundary condition. For irregular modes, 
a special case of boundary conditions at the cone apex is considered, when the MIT bag 
boundary condition is imposed at a finite radius, which is then taken to zero. We observe 
that the vacuum expectation values for both charge density and azimuthal current are peri- 
odic functions of the magnetic flux with the period equal to the flux quantum whereas the 
expectation value of the radial component vanishes. For both exterior and interior regions, 
the expectation values of the current are decomposed into boundary-free and boundary- 
induced parts. For a massless field the boundary-free part in the vacuum expectation value 
of the charge density vanishes, whereas the presence of the boundary induces nonzero charge 
density. Two integral representations are given for the boundary-free part in the case of a 
massive fermionic field for arbitrary values of the opening angle of the cone and magnetic 
flux. The behavior of the induced fermionic current is investigated in various asymptotic 
regions of the parameters. At distances from the boundary larger than the Compton wave- 
length of the fermion particle, the vacuum expectation values decay exponentially with the 
decay rate depending on the opening angle of the cone. We make a comparison with the 
results already known from the literature for some particular cases. 



PACS numbers: 03.70.-hk, 04.60.Kz, 11.27.-hd 



1 Introduction 

Topological defects are inevitably produced during symmetry breaking phase transitions and play 
an important role in many fields of physics. They appear in different condensed matter systems 
including superfluids, superconductors and liquid crystals. Moreover, symmetry breaking phase 
transitions have several cosmological consequences and, within the framework of grand unified 
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theories, various types of topological defects are predicted to be formed in the early universe [T]. 
They provide an important link between particle physics and cosmology. Among various types 
of topological defects, the cosmic strings are of special interest. They are candidates to produce 
a number of interesting physical effects, such as the generation of gravitational waves, gamma 
ray bursts and high-energy cosmic rays. Recently, cosmic strings attract a renewed interest 
partly because a variant of their formation mechanism is proposed in the framework of brane 
inflation [2]. 

In the simplest theoretical model describing the infinite straight cosmic string the spacetime 
is locally flat except on the string where it has a Dirac-delta shaped Riemann curvature tensor. 
Prom the point of view of quantum field theory, the corresponding non-trivial topology induces 
non-zero vacuum expectation values for several physical observables. Explicit calculations for 
the geometry of a single idealized cosmic string have been developed for different fields [3]- 
|24) . Moreover, vacuum polarization effects by higher-dimensional composite topological defects 
constituted by a cosmic string and global monopole are investigated in Refs. [25] for scalar and 
fermionic fields. The geometry of a cosmic string in background of de Sitter spacetime has been 
recently considered in |26j . 

Another type of vacuum polarization arises in the presence of boundaries. The imposed 
boundary conditions on quantum fields alter the zero-point fiuctuations spectrum and result in 
additional shifts in the vacuum expectation values of physical quantities. This is the well-known 
Casimir effect (for a review see [^27j). Note that the Casimir forces between material boundaries 
are presently attracting much experimental attention [28j. In Refs. j29]-[32j. both types of 
sources for the polarization of the vacuum were studied in the cases of scalar, electromagnetic 
and fermionic fields, namely, a cylindrical boundary and a cosmic string, assuming that the 
boundary is coaxial with the string. The case of a scalar field was considered in an arbitrary 
number of spacetime dimensions, whereas the problems for the electromagnetic and fermionic 
fields were studied in four dimensional spacetime. Continuing in this line of investigation, in the 
present paper we study the fermionic current induced by a magnetic fiux in a (2+l)-dimensional 
conical space with a circular boundary. 

As it is well known, field theoretical models in 2+1 dimensions exhibit a number of inter- 
esting effects, such as parity violation, flavour symmetry breaking, fractionalization of quantum 
numbers (see Refs. [33]-[39]). An important aspect is the possibility of giving a topological 
mass to the gauge bosons without breaking gauge invariance. Field theories in 2+1 dimensions 
provide simple models in particle physics and related theories also rise in the long-wavelength 
description of certain planar condensed matter systems, including models of high-temperature 
superconductivity. An interesting application of Dirac theory in 2+1 dimensions recently ap- 
peared in nanophysics. In a sheet of hexagons from the graphite structure, known as graphene, 
the long- wavelength description of the electronic states can be formulated in terms of the Dirac- 
like theory of massless spinors in (2+l)-dimensional spacetime with the Fermi velocity playing 
the role of speed of light (for a review see Ref. [30]). One-loop quantum effects induced by 
non-trivial topology of graphene made cylindrical and toroidal nanotubes have been recently 
considered in Refs. [H]. The vacuum polarization in graphene with a topological defect is 
investigated in Ref. |42j within the framework of long-wavelength continuum model. 

The interaction of a magnetic flux tube with a fermionic field gives rise to a number of 
interesting phenomena, such as the Aharonov-Bohm effect, parity anomalies, formation of a 
condensate and generation of exotic quantum numbers. For background Minkowski spacetime, 
the combined effects of the magnetic flux and boundaries on the vacuum energy have been 
studied in Refs. [431 144] . In the present paper we investigate the vacuum expectation value of 
the fermionic current induced by vortex conflguration of a gauge fleld in a (2+l)-dimensional 
conical space with a circular boundary. We assume that on the boundary the fermionic field 
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obeys MIT bag boundary condition. The induced fermionic current is among the most important 
quantities that characterize the properties of the quantum vacuum. Although the corresponding 
operator is local, due to the global nature of the vacuum, this quantity carries an important 
information about the global properties of the background spacetime. In addition to describing 
the physical structure of the quantum field at a given point, the current acts as the source in the 
Maxwell equations. It therefore plays an important role in modelling a self-consistent dynamics 
involving the electromagnetic field. 

From the point of view of the physics in the region outside the conical defect core, the 
geometry considered in the present paper can be viewed as a simplified model for the non-trivial 
core. This model presents a framework in which the influence of the finite core effects on physical 
processes in the vicinity of the conical defect can be investigated. In particular, it enables to 
specify conditions under which the idealized model with the core of zero thickness can be used. 
The corresponding results may shed light upon features of finite core effects in more realistic 
models, including those used for defects in crystals and superfluid helium. In addition, the 
problem considered here is of interest as an example with combined topological and boundary 
induced quantum effects, in which the vacuum characteristics can be found in closed analytic 
form. 

The organization of the paper is as follows. In the next section we consider the complete 
set of solutions to the Dirac equation in the region outside a circular boundary on which the 
field obeys MIT bag boundary condition. Shrinking the radius of the circle to zero we clarify 
the structure of the eigenspinors for the boundary-free geometry. These eigenspinors are used in 
Sect. [3] for the evaluation of the vacuum expectation value of the fermionic current density. Two 
integral representations are provided for the charge density and azimuthal component. In Sect, 
m we consider the vacuum expectation values in the region outside a circular boundary. They 
are decomposed into boundary-free and boundary-induced parts. Rapidly convergent integral 
representations for the latter are obtained. Similar investigation for the region inside a circular 
boundary is presented in Sect. El The main results are summarized in Sect. [H In Appendix 
[A] we derive two integral representations for the series involving the modified Bessel functions. 
These representations are used to obtain the fermionic current densities in the boundary-free 
geometry. In Appendix [Bl we compare the results of the present paper, in the special case of 
a magnetic flux in (2+l)-dimensional Minkowski spacetime, with those from the literature. In 
AppendixOwe show that the special mode does not contribute to the vacuum expectation value 
of the fermionic current in the region inside a circular boundary. 

2 Model and the eigenspinors in the exterior region 

In this paper we consider a two-component spinor fleld ip, propagating on a (2 + l)-dimensional 
background spacetime with a conical singularity described by the line-element 

ds'^ = Qf^udx^'dx" = df - dr^ - r'^d(p'^, (2.1) 

where r ^ 0, ^ (/> ^ and the points {r,4>) and (r, </> + (po) are to be identified. We are 
interested in the change of the vacuum expectation value (VEV) of the fermionic current induced 
by a magnetic flux in the presence of a circular boundary concentric with the apex of the cone. 
The dynamics of a massive spinor fleld is governed by the Dirac equation 

h^iV,, + ieA^)^ - = , = 9^ + r^, (2.2) 

where is the vector potential for the external electromagnetic fleld. In Eq. (j2.2p . = e|'^)7*-'^^ 
are the 2x2 Dirac matrices in polar coordinates and is the spin connection. The latter is 
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defined in terms of the flat space Dirac matrices, 7^"^^ , by the relation 

r, = {I'-'-h'^'^e'l^^e^.y,, , (2.3) 

where ; means the standard covariant derivative for vector fields. In the equations above, ^^^y 

a = 0,1,2, is the basis tetrad satisfying the relation 6(^)^(6)^"'' ~ 3^'^^ with t]""^ being the 
Minkowski spacetime metric tensor. We assume that the field obeys the MIT bag boundary 
condition on the circle with radius a: 

(l + m^7^)V|,^„ = , (2.4) 

where is the outward oriented normal (with respect to the region under consideration) to the 
boundary. In particular, from Eq. (12. 4p it follows that the normal component of the fermion 
current vanishes at the boundary, ni^tpj^t/j = 0, with ■0 = ip^'y^ being the Dirac adjoint and the 
dagger denotes Hermitian conjugation. In this section we consider the region r > a for which 

In (2+l)-dimensional spacetime there are two inequivalent irreducible representations of the 
Clifford algebra. In the first one we may choose the flat space Dirac matrices in the form 

^(0) = as, 7^^^ = iai, j^^^ = ia2, (2.5) 

with (T; being Pauli matrices. In the second representation the gamma matrices can be taken 
as 7^*^^ = — cr3, 7^^^ = —iai, 7^^^ = —ia2- In what follows we use the representation (|2.5p . The 
corresponding results for the second representation are obtained by changing the sign of the 
mass, m — )• —m. For the basis tetrads we use the representation below: 

^fi) = (0'Cos(#),-sin(#)/r) , (2.6) 
61*2) = (0,sin(#),cos(g(/))/r) , 

where the parameter q is related to the opening angle of the cone by the relation 

q = 27r/0o- (2.7) 



With this choice, for the Dirac matrices in the coordinate system given by the line element (|2.ip . 
we have the following representation 

Consequently, the Dirac equation takes the form 

7^ {d^ + ieA^) + ^^7^ +im]il) = 0, (2.9) 
2r J 

where the term with 7^ comes from the spin connection. 

In what follows we assume the magnetic field configuration corresponding to a magnetic flux 
located in the region r < a. This will be implemented by considering the vector potential in the 
exterior region, r > a, as follows 

A^ = {0,0,A), (2.10) 

In Eq. (j2.10p . ^2 = ^ is the covariant component of the vector potential in the coordinates 
{t, r, (j)). For the so called physical azimuthal component one has = —A/r. The quantity A is 
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related to the magnetic flux ^ by the formula A = —^/(pQ. Though the magnetic field strength, 
corresponding to ()2.10p . vanishes outside the flux, the non-trivial topology of the background 
spacetime leads to Aharonov-Bohm-like effects on physical observables. In particular, as it will 
be seen below, the VEV of the fermionic current depends on the fractional part of the ratio of 
^ by the quantum flux, 27r/e. 

Decomposing the spinor into upper and lower components, and cp-, respectively, from 
Eq. (12. 9j) we get the following equations 

(do ± im) v9± ± ie^^'i'f' \di + T -(^2 + ieA)] ip^ = 0. (2.11) 

L 2r r J 

From here we find the second-order differential equation for the separate components: 

.2 



dl - dl - \di - ^dl T 2i^d2 + % + m2 ) ^± = 0, (2.12) 



with the notations c± = (q — l)/2 it eA. 

For the positive energy solutions, the dependence on the time and angle coordinates is in 
the form e-*^*+*9"</' with > and n = 0,±1,±2, . . .. Now, from Eq. (|2T2]) . for the radial 
function we obtain the Bessel equation with the solution in the region r > a: 



where 7^0, 
and 

In Eq. ([2J3]1 . 



^+ = Z|,„|(7r)e^^'^^-^*, (2.13) 

^ = + m2, A„ = g(n + a + 1/2) - 1/2, (2.14) 
a = eA/q = -e<l>/2Tr. (2.15) 



Z\x„\{ir) = ciJ|A„|(7r) + C2l|A„|(7r), (2.16) 

with Ju{x) and Yjy{x) being the Bessel and Neumann functions. Note that in Eq. ()2.15p the 
parameter a is the magnetic flux measured in units of the flux quantum <&o = 27r/e. 

The lower component of the spinor, is found from Eq. (j2.1ip and for the positive energy 
eigenspinors we get 

V'S)(x) = e^^«<^-^*f JM^''^ , , V (2.17) 

where e^^ = 1 for A„ ^ and = — 1 for A„ < 0. From the boundary condition (j2.4p with 
= —6j^ we flnd 

^|An|(7«) + ;|^^|A„|+e,„(7a) = 0. (2.18) 
This condition relates the coefficients ci and C2 in the linear combination (I2.16P : 

C9 \ (7a) 

- = ^4—^- (2-19) 

Here and in what follows we use the notations (the notation with the upper sign is employed 
below) 

/W(z) = zf{z) + (±7^2^ ± ^ _ A„)/(z), /X = ma, (2.20) 
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for a given function f{z). 

Hence, outside a circular boundary the positive energy eigenspinors are presented in the form 



where 



_ „Jqn<p-iEt ( 5|A„|,|A„| (70,70 

V eA„i+;^5|A„MA„|+e,j7a,70 



(2.21) 



g,Ax,y) = Yy\x)Jp{y) - J^-\x)Yp{y). (2.22) 

Using the properties of the Bessel functions it can be seen that 

gv,u{x,y) = g-u,-u{x,y), g„,„+i{x,y) = -g-u-u-i{x,y). (2.23) 

Note that the spinor (j2.2ip is an eigenfunction of the operator J = —{i/q)d^ + CT3/2, with the 
eigenvalue j = n + 1/2, i.e., 

JiP\+\x) = jipi+Hx), j = n + 1/2. (2.24) 

The coefficient cq in Eq. ()2.2ip is determined from the orthonormalization condition for the 
eigenspinors: 

POO p(f>0 

J dr #rVft)^(x)^;t;,(x) = 5(7 - l')5nn' • (2.25) 

The integral over r is divergent when 7' = 7 and, hence, the main contribution comes from the 
upper limit of the integration. In this case, we can replace the Bessel and Neumann functions, 
having in the arguments the radial coordinate r, by the corresponding asymptotic expressions 
for large values of their argument. In this way, for the normalization coefficient we find, 

ter(7a)+lii:r(7a)l~\ (2.26) 



(l)o{E + m) L 



The negative energy eigenspinors are constructed in a similar way and they are given by the 
expression 

■j^(~)(2;) = CQe~^'^^'^~^^^^ I E+m fl'|A„|,|A„ 

with the same normalization coefficient defined by Eq. (|2.26|) . Note that the positive and 

negative energy eigenspinors are related by the charge conjugation which can be written as 

(— ) (+)* 

ip-yn = ciV'-yn , where the asterisk means complex conjugate. 

We can generalize the eigenspinors given above for a more general situation where the spinor 
field ip obeys quasiperiodic boundary condition along the azimuthal direction 

V^(t,r,,/. + ,/.o) = e2^^>^^(t,r,</)), (2.28) 

with a constant parameter x, |x| ^ 1/2. With this condition, the exponential factor in the 
expressions for the eigenspinors has the form e^*''^""'"^)'^^*-^* for the positive and negative energy 
modes (upper and lower signs respectively). The corresponding expressions for the eigenfunc- 
tions are obtained from those given above with the parameter a defined by 

a = x- e$/27r. (2.29) 

The same replacement generalizes the expressions for the VEVs of the fermionic current, given 
below, for the case of a field with periodicity condition (j2.28p . The property, that the VEVs 
depend on the phase x and on the magnetic flux in the combination (|2.29p . can also be seen by the 
gauge transformation = A'^ + d^A{x), ip{x) = ifj' {x)e~^^^^^\ with the function A(x) = A^x'^. 
The new function tp'{x) satisfies the Dirac equation with A'^ = and the quasiperiodicity 
condition similar to ()2.28p with the replacement X ^ x' = X ~ e^/^n. 
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3 Fermionic current in a boundary-free conical space 



Before considering the fermionic current in the region outside a circular boundary, in this section 
we study the case of a boundary-free conical space with an infinitesimally thin magnetic flux 
placed at the apex of the cone. The corresponding vector potential is given by Eq. (j2.10p for 
r > 0. As it is well known, the theory of von Neumann deficiency indices leads to a one-parameter 
family of allowed boundary conditions in the background of an Aharonov-Bohm gauge field [45j . 
In this paper, we consider a special case of boundary conditions at the cone apex, when the 
MIT bag boundary condition is imposed at a finite radius, which is then taken to zero (note 
that similar approach, with the Atiyah-Patodi-Singer type nonlocal boundary conditions, has 
been used in Refs. [44j for a magnetic flux in Minkowski spacetime). The VEVs of the fermionic 
current for other boundary conditions on the cone apex are evaluated in a way similar to that 
described below. The contribution of the regular modes is the same for all boundary conditions 
and the results differ by the parts related to the irregular modes. 

3.1 Eigenspinors 

In order to clarify the structure of the eigenspinors in a boundary-free conical space, we consider 
the limit a — )• for Eqs. (j2.2ip and (|2.27p . In this limit, using the asymptotic formulae for the 
Bessel functions for small values of the arguments, for the modes with j ^ —a we find 




E+m '^Pj+^j 



where 

^j=q\j + a\-ej/2. (3.2) 
Note that one has ej^j = Xn- In Eqs. (j3.ip and (|3.2p . we have defined 

' \ , (3.3) 

-1, J < -a ^ ' 

and the normalization coefficient is given by the expression 

(0)2 _ 7 E + m 



(3.4) 



In the case when a = N -\- 1/2, with being an integer number, the eigenspinors with 
j 7^ —a are still given by Eqs. (13. ip . The eigenspinors for the mode with j = —a, obtained from 
Eqs. (j2.2ip and (j2.27p in the limit a — )• 0, have the form (j3.ip with the replacements 

J^^.(z) ^ {E + m)Ji/2{z) -jYi/2{z), 
Jp^+e,{z) ^ (^ + m)J_i/2(z)-7y_i/2(z), (3.5) 

and ej = —1. The corresponding normalization coefficient is defined as c~q^ = {2E)~^ j/(j)Q. 
Taking into account the expressions for the cylinder functions with the orders ±1/2, the negative 
energy eigenspinors in this case are written as 

^0)7,—^^ \Trct)orEj \ e^?'^/^ cos(7r - 70) / 



7 



where 70 = arccos [y^{E — m) /2E] . Note that the eigenspinors obtained from in the hmits 
a —7- {N + 1/2)^ do not coincide with (j3.6p . For the hmit from below (above), a — )■ (A^ + l/2)~ 
(a —7- (N + 1/2)"'"), the eigenspinors are given by Eq. (|3.6|) with the replacement 70 — )• 7r/2 
(70 — ?• 0). Hence, for the bag boundary condition on the cone apex, the eigenspinors in the 
boundary-free geometry are discontinuous at points a = N + 1/2. Notice that, in the presence 
of the circular boundary, the eigenspinors in the region outside the boundary, given by Eqs. 
(I2.2ip and (]2.27p . are continuous. 

In general, the fermionic modes in background of the magnetic vortex are divided into two 
classes, regular and irregular (square integrable) ones. In the problem under consideration, for 
given q and a, the irregular mode corresponds to the value of j for which q\j + a| < 1/2. If we 
present the parameter a in the form 

a = ao + riQ, \ao\ < 1/2, (3.7) 

being hq an integer number, then the irregular mode is present if |ao| > (1 — l/q)/2. This mode 
corresponds to j = — no— sgn(ao)/2. Note that, in a conical space, under the condition 

|ao| ^ (1 - l/q)/2, (3.8) 

there are no square integrable irregular modes. As we have already mentioned, there is a one- 
parameter family of allowed boundary conditions for irregular modes, parametrized with the 
angle 6, ^ 9 < 27r (see Ref. [l5]). For |q!o| < 1/2, the boundary condition, used in deriving 
eigenspinors ()3.ip . corresponds to 9 = 3tt/2. If a is a half-integer, the irregular mode corresponds 
to j = —a and for the corresponding boundary condition one has 9 = 0. Note that in both cases 
there are no bound states. 



3.2 Vacuum expectation value of the fermionic current 

The VEV of the fermionic current, j^{x) = e^'j^ip, can be evaluated by using the mode sum 
formula 

poo 

(fix)) = eY, d^i>[-\xW^^-\x), (3.9) 
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where means the summation over j = ±1/2, ±3/2, .... In this section we consider this VEV 
for a conical space in the absence of boundaries. The corresponding quantities will be denoted 
by subscript 0. For the geometry under consideration, the eigenspinors are given by expressions 
(13. ip . Substituting them into Eq. p.9p one finds 



0-'(^))o = ^^Y.'^ r d^^^Jp^{^r)Jp^^,^{^r), (3.10) 

and the VEV of the radial component vanishes, (j^(x))o = 0. In deriving Eqs. (j3.10p . we have 
assumed that the parameter a is not a half-integer. When a is equal to a half-integer, the 
contribution of the mode with j = —a should be evaluated by using eigenspinors (j3.6p . The 
contribution for all other j is still given by Eqs. (|3.10p . As it will be shown below, both these 
contributions are separately zero and for half-integer values of a the renormalized VEV of the 
fermionic current vanishes. 
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In order to regularize expressions (|3.1U|) we introduce a cutoff function e~'^"' with the cutoff 
parameter s > 0. At the end of calculations the limit s — t- is taken. First let us consider the 
charge density. The corresponding regularized expectation value is presented in the form 



(/(^)>0,: 



reg 



eqm 
4tt 



E 



76 



-S7 



(3.11) 



Using the representation 



Jo 



(3.12) 



we change the order of integrations in the first term of the right-hand side in Eq. (j3.1ip and use 
the formula [36] 

d77e-^^'j|(7r) = -e^^'/^^I^ir^ /2s), (3.13) 

with Ip^z) being the modified Bessel function. The second term on the right of Eq. (|3.1ip is 
directly evaluated using Eq. (I3.13p . As a result, we get the following integral representation for 
the regularized charge density: 



0,reg 



eqme 



2(27r)3/2 



r^/2s ^-\l2p-m^r'^ I2z 

dz- 



Vr^ — 2zs 



+ 



eqe 



Svrs 



^[/^^.(rV2s)+/^^,+,^,(rV2.)] 



(3.14) 



The renormalization procedure for this expression is described below. 

Now we turn to the azimuthal component of the fermionic current. Using the relation 

zJ(}^+ej (z) = f3jJp^ (z) - ejzJ'p^ (z), (3.15) 

we write the corresponding regularized expression in the form 

f-oo e~'^'^^ Jg,{xr) 



rcg 



eg 

27rr2 



/■c 

^{e,Pj-rdr/2) / d7 7- 

■ Jo y'Y + 1^ 



(3.16) 



In a way similar to that we have used for the first term in the right-hand side of Eq. p. lip , the 
azimuthal current is presented in the form 



(/(^)>o, 



eqe 



E 



r2/2s 



dz 



^l/2^-m^r^/2z 
= — ( 

yr^ — 2zs 



[Ip,{z) - Ip^+^^iz)] . 



In deriving this representation we employed the relation 



(e,/3, - rdr/2) e'^'y Ip^{r'y) = ze"^ [l^^iz] - 1/3,+., (z)]^ 



(3.17) 



(3.18) 



for the modified Bessel function. 

The expressions of the regularized VEVs for both charge density and azimuthal current are 
expressed in terms of the series 

I{q,a,z)=J2l^A^)- (3-19) 

3 
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If we present the parameter a related to the magnetic flux as (j3.7|) . then Eq. (j3.19p is written 
in the form 

oo 

a, z) = ^ [^(n+ao + l/2)-l/2(^) + -fg(n-ao + l/2)+l/2(^;)] , (3-20) 
n=0 

which explicitly shows the independence of the series on ng- Note that for the second series 
appearing in the expressions for the VEVs of the fermionic current we have 

Y,lfi^^,^{z)=Z{q-a^,z). (3.21) 
i 

We conclude that the VEVs of the fermionic current depend on ao alone and, hence, these VEVs 
are periodic functions of a with period 1. 

When the parameter a is equal to a half-integer, that means |ao| = 1/2, the contribution 
to the VEVs from the modes with j 7^ —a is still given by Eqs. (I3.14p and ()3.17p . It is easily 
seen that for the case under consideration X^^y.Q, \l^.{x) — I^^.+e^. (x)] = 0. The contribution of 
the mode j = — a is evaluated by using eigenspinors (|3.6[) . A simple evaluation shows that this 
contribution vanishes as well. As regards the second term on the right-hand side of Eq. (|3.14p . 
below it will be shown that this term does not contribute to the renormalized VEV of the charge 
density. Hence, the renormalized VEVs for both charge density and azimuthal current vanish 
in the case when the parameter a is equal to a half- integer. Note that in the limit ao ~^ =tl/2, 
|ao| < 1/2, one has 

lim V (5X(g,5ao,2) =^\j2lT^ze~\ (3.22) 
and the expressions for the regularized VEVs are discontinuous at = ±1/2. 



3.3 Renormalized VEV in a special case 

Before further considering the fermionic current for the general case of the parameters charac- 
terizing the conical structure and the magnetic flux, we study a special case, which allows us 
to obtain simple expressions. It has been shown in [51 |9l [12] that when the parameter q is an 
integer number, the scalar Green function in four-dimensional cosmic string spacetime can be 
expressed as a sum of q images of the Minkowski spacetime function. Also, recently the image 
method was used in [22j to provide closed expressions for the massive scalar Green functions in 
a higher-dimensional cosmic string spacetime. The mathematical reason for the use of the image 
method in these applications is because the order of the modified Bessel functions that appear 
in the expressions for the VEVs becomes an integer number. As we have seen, for the fermionic 
case the order of the Bessel function depends, besides on the integer angular quantum number 
n = j — 1/2, also on the factor (g — l)/{2q) which comes from the spin connection. However, 
considering a charged fermionic field in the presence of a magnetic flux, an additional term will 
be present, the factor a. In the special case with q being an integer and 

a = l/2q-l/2, (3.23) 

the orders of the modified Bessel functions in Eqs. (j3.14p and (j3.17p become integer numbers: 
f3j = q\n\, j = n + 1/2. In this case the series over n is summarized explicitly by using the 
formula |46j 

f;' V(x) = _L|^e-cos(2.fc/9)^ (3 24) 

n=0 ^ k=0 

where the prime on the summation sign means that the term n = should be halved. 
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By making use of Eq. (j3.'24|) . for the regularized VEV of charge density we find 



p^pm'^s 1-} i'r^l2s -1/2 -m^r^/2z 

eme sr^ .->,,, s I , z e _^-2zsin^(Tvk/q) 



erne v— >, rr~/zs 
(/(x)>o.e, = (^Esin^(vr^/.)y^ dz'- 

+— ycos2(7r/t/(?)e"2(-V2.)sin2(.fe/,)^ (3^25) 

k=0 

In the hmit s — )• 0, the only divergent contribution to the right-hand side comes from the term 
k = 0. This term does not depend on the parameter q and is the same as in the Minkowski 
spacetime in the absence of the magnetic flux. Subtracting the k = term, then we take the 
limit s — 0. The integral is expressed in terms of the Macdonald function Ki/2{2mr sm^nk/q)) 
and for the renormalized charge density we find: 

9-1 

(/(^))o,ren = |^ sin(^A;/g)e-2"^^ ^-(^^Z") . (3.26) 

k=l 

Note that the second term on the right-hand side of Eq. (|3.14|) does not contribute to the 
renormalized VEV. Expression ()3.26p coincides with the result of Ref. [24j obtained by using 
the Green function approach. 

In a similar way, using (j3.24p . for the regularized VEV of the azimuthal current we get the 
formula 

m^s fr^/2s 1/2 ~m'^r'^/2z 

(/(x))o,reg = ^^^sinHnk/q) / dz ' e-^-'° (3.27) 

vrr^ V 2vr ^ Jo Vr^ - 2zs 

This expression is finite in the limit s — )• and for the renormalized VEV one finds 

ovrr-^ ^— ^ sm[TTk/q) 

After the coordinate transformation (p' = qcf), this expression coincides with the result given in 
Ref. [24] ■ Note that for both charge density and azimuthal current one has (j'^(a;))o.rcn/e ^ 0. 
As expected the renormalized current densities decay exponentially at distances larger than 
the Compton wavelength of the fermionic particle. In Fig. [T] the VEVs of charge density and 
azimuthal current are plotted versus mr for different values of q. The corresponding values of 
the parameter a are found from Eq. (|3.23p . 



3.4 General case 

Now we turn to the general case for the parameters q and a. As it follows from Eqs. (j3.14p and 
p.lTh . we need to evaluate the integrals 

rr^/2s ^± 1/2 g-m^r^ /2a; 

/ dx e~''[I{q,ao,x) -I{q,-ao,x)], (3.29) 

JO V?' — 2xs 

in the limit s — )• 0, with the function I{q, aQ,x) defined by Eq. (]3.19p . Two alternative integral 
representations for this series are given in Appendix El In order to provide an integral represen- 
tation for the VEV of fermionic current, we first consider the representation ()A.6p . In the limit 
s — ?■ 0, the only divergent contributions to the integrals in Eq. (I3.29P with separate functions 
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mr mr 



Figure 1: VEVs of the charge density (left plot) and the azimuthal current (right plot) in a 
boundary-free conical space, as functions of the parameter mr for the special case of integer 
values of q with the magnetic flux defined by Eq. (|3.23p . 



X(g, ibao,x) come from the first term in the right-hand side of Eq. ()A.6p . This term does not 
depend on oq and, consequently, it is cancelled in the evaluation of the integral in Eq. (|3.29|) . 

Hence, we see that the regularized expression (j3.17p for the VEV of azimuthal current is 
finite in the limit s — >• 0. For the corresponding renormalized quantity we find 

3 /"OO 

(j'(x))o,ren = dz Z^l\-^''' '''^-^ J] bX{qM^,z). (3.30) 



5=±1 



where the function X((7, ao,^;) is given by the integral representation (lA.6p . The VEV of the 
azimuthal current is a periodical function of the magnetic flux with a period equal to magnetic 
flux quantum $o- It is an odd function of the parameter oq defined by Eq. (13. 7p . From Eq. (IA.6P 
it is seen that the integrand in Eq. (j3.30p decays exponentially for large values of z. Substituting 
the integral representation (|A.6P into Eq. (j3.30p and changing the order of integrations in the 
part with the last term of Eq. ()A.6p , the integral over z is expressed in terms of the Macdonald 
function K^i2{y). 

As a result, for the renormalized VEV of the azimuthal component we find the following 
expression 

2/ e r (— 1)' sin(27r/ao) 1 + 2mr sin(7r//g) 

_q_ r ^ E5=±i^fiQ^^(^0^y) l + 2mrcosh(j//2) -| 
47r io cosh(g'y) - cos(g7r) cosh^(y/2)e2''"''cosh(j//2) J ' 

where p is an integer defined by 2p < g < 2p -|- 2 and the function f{q,ao,y) is given by Eq. 
([ATSl) . In the case q = 2p, the additional term 

- (-1)^/2 sin(7rgao)e-2™'^(l/2 + mr), (3.32) 

should be added to the expression in the figure braces of Eq. (j3.3ip . For 1 ^ q < 2, only the 
integral term remains. The difference of the functions f{q,aQ,y), appearing in the integrand, 
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can also be written in the form 



6 f{q,6ao,y) = 2 coshiy/2) ^ (5 cos [gvr (1/2 - Jao)] cosh [g (1/2 + (5ao) y] . (3.33) 
5=±l 5=±l 

Note that for q = 2p the integrand in the last term of Eq. (j3.3ip is finite at y = 0. 

In the massless limit the expression in the figure braces of Eq. (I3.3ip does not depend on 
the radial coordinate r and the renormalized VEV of the azimuthal cm'rent behaves as 1/r'^. 
For a massive field, at distances larger than the Compton wavelength of the spinor particle, 
mr ^ 1, the VEV of the azimuthal current is suppressed by the factor e~^™^' for 1 ^ q < 2 and 
by the factor g-Sm-'" sin(7r/g) £qj. q ^ 2. In the limit mr <^ 1, the leading term in the corresponding 
asymptotic expansion coincides with the VEV for a massless field. In Fig. [2] we plot the VEV of 
the azimuthal current for a massless fermionic field as a function of the magnetic flux for several 
values of the parameter q (numbers near the curves). In the limit ag — ±1/2, |ao| < 1/2, the 
VEV of the azimuthal current is obtained using relation (|3.22p . From Eq. (|3.30p one finds: 

lim {j\x))o,ren = T^^Ki{2mr), (3.34) 
where K^(z) is the Macdonald function. 



4 
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Figure 2: The VEV of the azimuthal current in a boundary- free conical space as a function of 
the magnetic flux for a massless fermionic field. 

Now we turn to the VEV of the charge density. In the corresponding regularized expression 
(I3.14p . the first term on the right-hand side with the integral is finite in the limit s — t- 0. The 
second term is written in the form 

^S^Z^ V I{q, 5aoy/2s). (3.35) 

Taking into account Eq. ()A.6p . we see that the only nonzero contribution to this term comes 
from the first term in the right-hand side of Eq. (|A.6p . This contribution diverges and does not 
depend on the parameters q and a^. The divergence is the same as in the Minkowski spacetime 
in the absence of the magnetic flux and is subtracted in the renormalization procedure. As a 
result, for the renormalized VEV of the charge density we get the formula 

(J°(^))0,re„ = TTTTtS^ r .-^/2e~"''-V2.-. ^ ^^(^^ ^^^^ (3.35) 

2(27r)'^/^r 
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This VEV is a periodical function of the magnetic flux with a period equal to the magnetic flux 
quantum. For 2p < q < 2p -\- 2, with p being an integer, using the integral representation ()A.6p . 
the renormalized VEV is presented in the form 



1=1 

q r ^y ^S=±i SfiQ, Sao,y) g-^'nr coshfa/2) 
47r Jq cosh(gy) — cos((77r) cosh(y/2) J 
When q = 2p, the additional term 

- (-1)9/2 sin(7rgao)e-2"'V2, (3 3g) 

should be added to the expression in the figure braces on the right-hand side of Eq. ()3.31|) . As 
in the case of the azimuthal current, the renormalized VEV (j3.36p is an odd function of the 
parameter oq. Note that in a boundary- free conical space the charge density for a massless 
field vanishes at points outside the magnetic flux. For a massive field, at distances larger than 
the Compton wavelength, mr ^ 1, the renormalized charge density is suppressed by the factor 
g-2mr for 1 ^ g < 2 and by the factor e-^'^'-'^i"^'^/'') for q^2. 

Though the charge density given by Eq. (13.360 diverges at r = 0, this divergence is integrable 
and the total fermionic charge 

Q = (t)Q i drr(j°(a;))o.ren = 7 / dxe'"^ V 5X{q,5ao,x), (3.39) 
Jo 4 Jo ^^^^ 

is finite. In the form given by Eq. (|3.39|) . the integrals with 5 = 1 and 5 = —1 diverge 
separately and we cannot change the order of the integration an summation over 5. In order to 
overcome this difficulty, we write the integral as dx • • • = limg_^|+ dx e~*^' • • • . With 
this representation, evaluating the integrals with separate 6 for s > 1 and taking the limit s — )• 1, 
one finds 

Q = -eao/2. (3.40) 

This result for a conical space was previously obtained in Ref. [14]. As we see, the total charge 
does not depend on the angle deficit of the conical space. This property is a consequence of the 
fact that the total charge is a topologically invariant quantity depending only on the net flux. 

In the absence of the angle deficit one has q = 1 and the expressions for the VEVs of the 
fermionic current are simplified to 

_ emsin(7rao) ,cosh{2aoz) nmrcoshz 



= ^^^^./o coshz ' 



(j^(x))oren = _ ^^Wvrao) /-^^ cosh(2ao.) + ^mr cosh z) e-^— (3.41) 
47r^r'^ Jo cosh'^z 

Alternative expressions for the VEVs of the charge density and azimuthal current in (2+1)- 
dimensional Minkowski spacetime in the presence of a magnetic flux were given in Ref. [35j (see 
also |38j for the general case of a one-parameter family of boundary conditions at the origin). 
We compare these expressions with the formulae obtained in the present paper in Appendix [Bl 
In Fig. [31 the VEVs of the charge density (left panel) and azimuthal current (right panel) 
are plotted as functions of the magnetic flux for a massive fermionic field in a conical space with 
(po = TT. In the limit oq — >• ±1/2, |ao| < 1/2, for the azimuthal current one has Eq. (I3.34p . For 
the charge density the limiting values are given by 

hm {f{x))o,en = T^Ko{2mr). (3.42) 
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Both charge density and azimuthal current exhibit the jump structure at half- integer values 
for the ratio of the magneitc flux to the flux quantum (for a similar structure of the persistent 
currents in carbon nanotube based rings see, for example, Refs. [47J). 




Figure 3: The VEVs of the charge density (left panel) and the azimuthal current (right panel) 
as functions of the magnetic flux for a massive fermionic field in a boundary-free conical space 
with q = 2. 

Alternative expressions for the VEVs of the fermionic current are obtained by using the 
integral representation (jA.lSP for the functions X[q,^aQ, z). We start with the charge density. 
The corresponding regularized expression is given by Eq. ()3.14p . As we have already noticed, 
the first term on the right of this formula is finite in the limit s — )• 0. After the application 
of (jATTHI) to the series in the second term, we see that the parts corresponding to the first and 
last terms in the right-hand side of ()A.15P vanish in the limit s — >• due to the exponential 
decay of the Macdonald functions. The only term which survives in the limit s — )• is the part 
corresponding to the second term in the right-hand side of Eq. (jA.lSP . In the expression for the 
fermionic current this term is multiplied by q and hence it does not depend on the angle deficit 
and on the magnetic flux. So, this term is the same as in the case of Minkowski spacetime in the 
absence of the magnetic flux and is subtracted in the renormalization procedure. As a result, 
for the renormalized charge density one finds the expression below: 

sgn(ao)gS((?(|ao|-l/2) + l/2,z) + 2 / dyKiy{z)g{q,ao,y) , (3.43) 

JO ^ 



where 



B{y,z) = { . , ° ' ^ , , (3.44) 



and we have defined the function 



9iq,ao,y) = X] 
<5=±1 



5 sinh(y7r) 



p27r(t/+i|g5ao-l/2l)/g _^ ]_ 



(3.45) 
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For the VEV of the azimuthal current, m a similar way, we get the representation 



4er-3 f 



oo 



,^r^ /2z—z 



,ren — 



(27r)5/2 I 



oo 



X sgn{ao)qB{q{\ao\ - 1/2) + 1/2, z) + 2 



dy Kiy{z)g{q,ao,y) 



(3.46) 



Note that under the condition (j3.8p there are no square integrable irregular modes and, in this 
case, the first terms in the square brackets of Eqs. (|3.43|) and (|3.46p vanish. 

In the special case q = 1 we see that g{l,ao,y) = 0, and for the VEVs of the fermionic 
current we obtain the formulae 



In the limit oq — )■ ±1/2 we recover results (lOil) and (13:421) . In Appendix El we show the 
equivalence of these expressions to the ones previously given in the literature for the fermionic 
densities induced by a magnetic flux in (2+l)-dimensional Minkowski spacetime. 

In the discussion above we used the irreducible representation of the Clifford algebra corre- 
sponding to Eq. ()2.5p . For the second representation the renormalized VEV of the azimuthal 
current is given by the same expressions, whereas the expressions for the VEV of the renormal- 
ized charge density change the sign. Consequently, the total induced charge (j3.40p changes the 
sign as well. 

4 Induced fermionic current in the exterior region 

Now we turn to the investigation of the induced fermionic current in the presence of a circular 
boundary at r = a with boundary condition (j2.4p . The corresponding VEV is evaluated using 
the mode sum formula ()3.9p with the eigenspinors given by Eq. (I2.27P . For the further discussion 
it is convenient to write these eigenspinors in an equivalent form by using the properties ()2.23p : 



where cq is given by Eq. ()2.26p with the replacement | A„ | — )• /3j . In Eq. (j4.ip , ej is defined by 
(|3.3p for j ^ —a and ej = —1 for j = —a. We could also obtain the representation ()4.ip . by 
taking instead of Eq. ()2.16p the linear combination of the functions J/Sji'yr) and Yg^. (7r). Note 
that the corresponding barred notation in expression (12.221) may also be written in the form 




(3.47) 




(4.1) 




ejzF^^+,^ (z) ± (^^2 + ^2 + ^)^^^ 



(4.2) 



with F = J,Y and fi = ma. 
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As in the boundary-free case, the VEV of the radial component vanishes and for the charge 
density and the azimuthal current we have the expressions 



(i°(^)> 



eg 



j{E- m)gj^^p^^^^ (7a, 7r) + {E + m)gj. p^ (7a, jr) 



E 



^(-)2, 







E 4-)^(7a) + F^-)^(7a) 



(4.3) 



Here, as before, the summation goes over j = ±1/2, ±3/2, . . .. Equivalent forms are obtained 
using the eigenspinors ()2.27p . From Eq. (j4.3p it follows that the fermionic current is a periodic 
function of a with the period equal to 1. We assume that a cutoff function is introduced without 
explicitly writing it. The specific form of this function is not important for the discussion below. 
In the presence of the boundary, the VEV of the fermionic current can be decomposed as 



{f{x)) = {f{x))o + ifix))^, 



(4.4) 



where {j'^{x))\y is the part induced by the boundary. In order to extract the latter explicitly, we 
subtract from Eq. (j4.3p the VEVs when the boundary is absent. If the ratio of the magnetic 
flux to the flux quantum is not a half- integer, the boundary- free parts are given by Eq. (I3.10p . 
In this case, for the evaluation of the difference, in the expression of the charge density we use 
the identity 

.1^ [ 'ID = — ^ 



(-)2. 



Jl{y) 



Jj-\x) 



H'~:''{x) 



1=1,2 ^^/3, 



(4.5) 



with A = /3j,/3j ± Cj, and with the Hankel functions He, {x). The expression for the boundary- 
induced part in the VEV of the charge density takes the form 



(i°(^)>b 



eq_ 

'Svr 



EE 

j 1=1,2 

M)2 



7 "^/3 



4:^(7«) 



Efji-'i) 



{E - m)Hi;iJjr) + {E + m)H'^>\-fr) 



(7a) 

(02, 



(4.6) 



Now, in the complex plane 7, we rotate the integration contour by the angle it/2 for the term 
with I = 1 and by the angle — 7r/2 for the term with / = 2. The integrals over the segments 
(0, im) and (0, —im) cancel each other and, introducing the modified Bessel functions, we get 
the following expression 



(/(^)>i 



eq_ 
'2^2 



E 



dz z 



X < m 



Re 



I^-\za)/Ky/{za) 



(')< 



Im 



4)(za)/4^,)(za)] }, 



where we use the notation (the notation with the upper sign is used in the next section) 

f(±)(. 



zF'{z) + ( ±^ ± iVz'^ - /^^ 



e,f3AF{z). 



(4.7) 



(4.8) 
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The ratio in the integrand of Eq. (|4.7|) can also be written in the form 



(4.9) 



with the notation 



Wl,%^+,^ (x) = X [Ip^ {x)Kp^ (x) - Ip^^,^ {x)K^,+e, (x)] 

±fi [1^^,+,^, {x)Kp^ (x) - Ip^ {x)Kp^+,^ {x)] . (4.10) 

The real and imaginary parts appearing in Eq. (j4.7p are easily obtained from Eq. (j4.9p . Note 
that under the change a — )■ —a, j — —j, we have /3j — )• + e^, /3j + ej — )• /3j. From here it 
follows that the real/imaginary part in Eq. ()4.9p is an odd/even function under this change. 
Now, from Eq. (j4.7p we see that the boundary-induced part in the VEV is an odd function of 
a. When a is a half-integer, in the term of Eq. (j4.7p with j = —a the orders of the modified 
Bessel functions are equal to ±1/2. Using the expressions of these functions in terms of the 
elementary functions, we can see that the corresponding integral vanishes. For the terms with 
j 7^ —a, the contributions of j < —a and j > —a to the right-hand side of Eq. (j4.7p cancel each 
other. Hence, if a is a half-integer the boundary-induced part in the VEV of the charge density 
vanishes. Recall that the same is the case for the boundary-free part. 
In a similar way, using the identity 

g,^.(x,,)g,,.^ (X,,) ^ ^ _ 1 y ^k!^Hl\y)HS\ (y)^ (4-11) 

for the boundary-induced part in the VEV of the azimuthal current we find 

v2 



xKp^{zr)Kis^+,^{zr)Re \lj;-\za)/Kj,-\za)\ . (4.12) 



As in the case of charge density, this part is a periodical function of the magnetic flux with a 
period equal the flux quantum. 

If we present the ratio of the magnetic flux to the flux quantum in the form (|3.7p . then the 
boundary induced VEVs are functions of oq alone. They are odd functions of this parameter. 
In the limit ao ~^ |ao| < 1/2, the only nonzero contribution to (i^(a;))b comes from the 

term with j = ±1/2 and one has the limiting values 

lim (/(x))b = ±^Koi2mr), 

lim ifix))^ = ±^K,{2mr). (4.13) 

Now comparing with Eqs. ()3.34p and (j3.42p . we see that the limiting value of the total current 
density (|4.4p is zero and the latter is continuous at oq = ±1/2. This result was expected due 
to the continuity of the exterior eigenspinors as functions of the parameter oq. Comparing with 
the results of the previous section, we see that the limiting transitions a — t- and |ao| — 1/2 do 
not commute. 
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For a massless field the expressions for the boundary-induced parts in the VEVs take the 
form 



(J°(^))b 
(j'(^))b 



eq 



27r2a2 



E 



°° , KliMr/a)-K' {zr/a) 
dz 







eq 



~ Kp^{zr/a)Kp^+,^{zr/a) ^^) 







(4.14) 



with the notation defined by Eq. (j4.10p . We would like to point out that the boundary-induced 
charge density does not vanish for a massless filed. The corresponding boundary-free part 
vanishes and, hence, {j^{x)) = {j^{x)))^. 

Now we turn to the investigation of the boundary-induced part in the VEV of fermionic 
current in the asymptotic regions of the parameters. In the limit a — t- 0, for fixed values of r, 
by taking into account that 



- a 



2 — mP' 



k[ \za 

Pi 



r2(g|i + a| + l/2) 

to the leading order, from Eqs. ()4.7p and (j4.12p . we have 
eq sgn(ao)(a/2r)2'Ja 



za 



I2f<i\j 



+Ol\-l 



(4.15) 



(/(^)>b ^ 

0-'(^)>b ■ 

with the notation 



7r2 r2r2(g„ + l/2) 



2qa 



dz 



\/ z^ — vn?r'^ 



(2mV2 - z2) Kl^_^i^{z) + z^Kl^^^i^{^z) 



(4.16) 



2eqm sgn(Qo)(a/2r)2'?° 
7r2r2 r2(g^ + l/2) 



dz 



z 



\/ z'^ — m'^r'^ 



g„ = g(l/2 - |ao|). 



(4.17) 



For a massless field the asymptotic behavior for the charge density is directly obtained from Eq. 
(j4.16p . The integrals involving the Macdonald function are evaluated in terms of the gamma 
function and one finds 



0/ XV egsgn(ao) / a \ 29a qaT{2qa + l/2)r(ga + 1) 
OW)b~ \2r) (2g, + l)r3(g„ + 1/2) 



(4.18) 



For the azimuthal component the leading term in Eq. (j4.16p vanishes. The corresponding 
asymptotic behavior is directly found from Eq. (j4.14p . The leading term is given by 



0-'(^))b 



2eq sgn(ao) / a \ 2ga+i r(2g„ + 3/2)r(g„ + 1) 



7rr3 (2(?„)2 - 1 \2r) {2q^ + l)r3(g„ + 1/2) ' 

for qa > 1/2, and by the expression 

egsgn(ao) / a \ 4ga r(l/2 - g«) 



0"(^)>b 



-r(2g« + i/2)r(3(7a + 1) 



22'?c.+i7rr3 \2rJ T^{l/2 + qa) 
in the case qa < 1/2. For for qa = 1/2, the leading terms behaves as a^ln(a) 



(4.19) 



(4.20) 
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At large distances from the boundary, for a massive field, under the condition mr ^ 1, 
the dominant contribution to the integrals come from the region near the lower limit of the 
integration and to the leading order we find 



(J°(^)>b 
(j'(^)>b 



eqmPe 
4-y/7r(mr)^/2 



^Re [/(;)(ma)/4;)( 
j 

5; Re [4-)(ma)/4-)( 



(4.21) 



As expected, we have an exponential suppression of the boundary-induced VEVs. For a massless 
field, the asymptotics at large distances are given by Eqs. (|4.18p - ()4.20p . In Fig. HI we plot the 
VEVs of the charge density (left panel) and azimuthal current (right panel) for a massless 
fermionic field as functions of the magnetic flux. The graphs are plotted for r/a = 1.5 and for 
several values of the parameter q. As we have already mentioned, in the exterior region the 
total VEVs of the charge density and azimuthal current vanish in the limt |ao| ~^ 0. Note that 
for a massless field the boundary-free part in the VEV of the charge density vanishes and the 
non-zero charge density on the left plot is induced by the circular boundary. 
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Figure 4: The VEVs of the charge density (left panel) and azimuthal current (right panel) for a 
massless fermionic field outside a circular boundary. 



5 Fermionic current inside a circular boundary 

In this section we consider the region inside a circular boundary with radius a, r < a, on which 
the fermionic field obeys the boundary condition ()2.4p with = —6j^. The boundary condition 
at the cone apex for the irregular mode is the same as that we have used in Section [3] for the 
boundary-free conical geometry. The eigenspinors in this region have the form 



= ^oe-«<^+^^*f 1, (5.1) 



\-m 

iqrp/2 



^3^ 

E+m '^Pj+^j 



e 
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with the same notations as in Section [31 From the boundary condition at r = a we find that 
the eigenvalues of 7 are solutions of the equation 



(5.2) 



Note that this equation may also be written in the form j'^^^ (7a) = with the barred notation 
defined by Eq. (I2.20p . For a given 13 Eq. (I5.2p has an infinite number of solutions which we 
denote by 7a = 7^^.,;, / = 1, 2, . . .. 

The normalization coefficient in Eq. (15. ip is determined from the condition 



dr 



00 



Using the standard integral for the square of the Bessel function (see, for example, 
finds 



2(y2 + _ + 1) ^y2 + ^2 + 



, y = 7/3,,/, 



(5.3) 
I, one 

(5.4) 



where, as before, /U = ma. For the further convenience we will write this expression in the form 



9^0 



yTp, jy) /i + _^ ^2 



with the notation 



y2 + (/, - e^./3j.) U + 



(5.5) 



(5.6) 



Substituting the eigenspinors (15. ip into the mode sum formula 



(/(^)> = e^^v5;7^(x)7>;7(3 
j 1=1 

for the VEVs of separate components of fermionic current we have 



(5.7) 



eg 



j 1=1 



+ l]jl{yr/a) 



1 J; 



27r a^r 



JlSj{yr/a)Jfs^+,^{yr/a), 



2 



(yr/a) 
(5. 



with y = 7/3j,i, and the radial component vanishes, (j^(x)) = 0. As before, we assume the 
presence of a cutoff function without explicitly writing it. 

As the explicit form for 'jpj^i is not known, Eqs. (15. 8p are not convenient for the direct 
evaluation of the VEVs. In addition, the separate terms in the mode sum are highly oscillatory 
for large values of the quantum numbers. In order to find a convenient integral representation, 
we apply to the series over / the summation formula (see 



^ poo poo 

y]/(7/3,,;)^/3(7/3i,;) = / dxf{x)-- dx 

Jo ^ JO 

K^-^^ix) K^'^^*(x) 
-ft-"7(xe"/^) ^\ V e^^"/(xe-"/^)^ 



(5.9) 
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the asterisk meaning complex conjugate. Here the notation F^'^^x) for a given function F{x) 
is defined by Eq. (j4.8p for x ^ ^ and by the relation 

= xF'{x) + {fi+ y^n^ - x^ - ejf3j)F{x), (5.10) 

for X < iJ,. Note that in the latter case F^^^*{x) = F^'^^x). The term in the VEVs corresponding 
to the first integral in the right-hand side of Eq. (j5.9p coincides with the VEV of the fermionic 
current for the situation where the boundary is absent. 

As a result, the VEV of the fermionic current is presented in the decomposed form (j4.4p . 
For the function f{x) corresponding to Eq. (15. Sp . in the second term on the right-hand side of 
Eq. (j5.9p . the part of the integral over the region (0,^) vanishes. Consequently, the boundary- 
induced contribution for the charge density in the region inside the circle is given by the expres- 
sion 

oo 



o-V))b = -|^E/J^-- 



f /i (zr) + /i , ^(zr) 
x{m ' MK^^\za)/I^;^\za)] 

-Vlizr) - ll^^^{zr)]lr^[K^^\za)/I^^\za)]]. (5.11) 

Similarly, for the boundary-induced part in the azimuthal component we find 

For points away from the circular boundary, the boundary-induced contributions (jS.lip and 
(I5.12P are finite and the renormalization is reduced to that for the boundary-free geometry. 
These contributions are periodic functions of the parameter a with the period equal to 1. So, 
if we present this parameter in the form (13. 7p with no being an integer, then the VEVs depend 
on ao alone and they are odd functions of this parameter. Similar to Eq. ()4.9p . the ratio of the 
combinations of the modified Bessel functions in Eq. ()5.1ip is presented in the form 



4+) {x) ^1% (^) + ^1,+., (^)] + 2m^/3, (x)//3,+., (x) - 



(5.13) 



with the notation defined by Eq. (j4.10p . 

When the parameter a is a half-integer the contributions of the modes with j ^ —a to the 
boundary-induced VEVs inside the circle are still given by expressions (jS.lip and (j5.12p . It 
can be easily seen that the contributions of the modes with j < —a and j > —a cancel each 
other. The contribution of the mode with j = —a should be considered separately. In Appendix 
Owe show that this contribution vanishes as well. Therefore, we conclude that for a being a 
half-integer, the boundary-induced part in the VEV of the fermionic current vanishes. 

For a massless field the formulae of the boundary-induced parts are simplified to 

r-oo Il{zr/a)-ll.(zr/a) 



(fix)), - Td.^^^^^^^^^— ^^^^ 
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Note that for a massless field VF^^^.^^.(z) = ^.^^.(-z). In the hmit ao — )• ±1/2, |ao| < 1/2, 
the only nonzero contributions to Eqs. (|5.1ip and (|5.12|) come from the term with j = ^1/2 
and, by making use of Eqs. ()3.34p . for the total VEVs we find 

eqm \ _s_ ^9 f°° j amzcosh(2zr) — (z + m)e^^" 



hm {fix)) = ^^Ko{2mr)±^ r dz 



eg^m ega /""^ sinh(2zr) 



For a massless field they reduce to the expressions: 



eq 



lim (i^(x)) = =F , 
lim (i^(x)) - - 



7rr/2a 
1+ ' 



sin(7rr/2a) 



(5.16) 



Note that the limiting values are linear functions of the parameter q. 

The general expressions for the VEVs are simplified in asymptotic regions of the parame- 
ters. First we consider large values of the circle radius. For the modified Bessel functions in 
the integrands of Eqs. ()5.1ip and ()5.12p . with za in their arguments, we use the asymptotic 
expansions for large values of the argument. By taking into account that for a massive field the 
dominant contribution into the integrals comes from the integration range near the lower limit, 
to the leading order we find 

2 —2ma 

ifix))^ « _^«^L^^^(2e^./3. + l)I (^,)/ (^,). (5.17) 

In this limit, for a fixed value of the radial coordinate, the boundary-induced VEVs decay 
exponentially. For a massless field, assuming r/a ^ 1, we expand the modified Bessel function 
in the numerators of integrands in Eq. (j5.14p in powers of r/a. The dominant contribution 
comes from the term j = 1/2 for oq < and from the term j = —1/2 for ao > 0. To the leading 
order we have 

/■Of A\ eg sgn(ao)(r/2a)2''"-i ^'^""^ 

vr^a^ {2q^ + l)T\qo, + 1/2) Jq Iq^+i/2^z) + ^^-1/2^^) 

where q^ is defined in Eq. ()4.17p . As it is seen, for a massless field the decay of the VEVs is as 
power-law. 
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For points near the apex of the cone, r — )• 0, we have the following leading terms 



(J°(^))b 



2qa — l roo 



eg sgn(ao)(r/2a) 
27r2a2 r2(g„ + l/2 



dz 



qa-1/2 



z)]+2nlg^-l/2{z)lq^ + l/2i 



(j'(x))b 



eg sgn(ao)(r/2a) 
'7r2a3 (2q„ + l)r2(g, + 1/2) 

qQ-l/2,qQ + l/2^ 



2; 



dz 



z 



2qc^2 



^[^l-l/2(^) + ^i + l/2(^)] + Wg.-l/2(^)^„ + l/2(^ 



(5.19) 



(5.20) 



For a massless field these expressions reduce to Eq. (|5.18|) . From here it follows that in the 
limit r — )• the boundary- induced part vanishes when |ao| < 1/2 — l/(2g) and diverges for 
|ao| > 1/2 — l/(2g). Notice that in the former case the irregular mode is absent and the 
divergence in the latter case comes from the irregular mode. This divergence is integrable. In 
the case |q;o| = 1/2 — 1/(2(7), corresponding to ga = 1/2, the boundary-induced VEV tends 
to a finite limiting value. In particular, for the magnetic vortex in the background Minkowski 
spacetime, the boundary-induced contribution diverges as r~2l"ol_ jjj Yig. [5l the VEVs of 
the charge density (left panel) and azimuthal current (right panel) are plotted for a massless 
fermionic field inside a circular boundary as functions of the magnetic flux. The graphs are 
plotted for r/a = 0.5 and for several values of the opening angle for the conical space. For a 
massless field the boundary-free part in the VEV of the charge density vanishes and the charge 
density on the left plot is induced by the boundary. 
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Figure 5: The same as in figure H] for the region inside a circular boundary. 



6 Summary and conclusions 

In this paper we have investigated the VEV of the fermionic current induced by a magnetic 
fiux string in a (2-|-l)-dimensional conical spacetime with a circular boundary. The case of 
massive fermionic field is considered with MIT bag boundary condition on the circle. In (2-1-1)- 
dimensional spacetime there are two inequivalent irreducible representations of the Dirac matri- 
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ces. We have used the representation (|2.5p . The corresponding results for the second represen- 
tation are obtained by changing m — )■ —m. Under this change, the boundary-free contribution 
to the VEV of the azimuthal current is not changed, whereas the VEV of the charge density 
changes the sign. For the evaluation of the expectation values we have employed the direct mode 
summation method. The corresponding positive and negative energy eigenspinors in the region 
outside the circular boundary are given by Eqs. (|2.2ip and (j2.27p . 

The VEV of the fermionic current in the boundary-free conical space is investigated in Sect. 
[3l For this geometry, under the condition (|3.8p . there are no square integrable irregular modes. 
In the case |qo| > (1 — l/q)/2, the theory of von Neumann deficiency indices leads to a one- 
parameter family of allowed boundary conditions at the origin. Here we consider a special 
boundary condition that arises when imposing bag boundary condition at a finite radius, which 
is then shrunk to zero. The VEVs of the fermionic current for other boundary conditions on the 
cone apex are evaluated in a similar way. The contribution of the regular modes is the same for 
all boundary conditions and the formulae differ by the parts related to the irregular modes. For 
the boundary condition under consideration, the eigenspinors for the boundary-free geometry 
are obtained from the corresponding functions in the region outside a circular boundary with 
radius a, taking the limit a — )• 0. They are presented by Eq. (j3.ip . When the magnetic flux, 
measured in units of the flux quantum, is a half-integer, there is a special mode corresponding 
to the angular momentum j = —a, with the negative-energy eigenspinor given by Eq. (j3.6p . 

In the boundary-free geometry, the regularized expressions of the VEVs are given by Eqs. 
(j3.14p and (j3.17p for the charge density and the azimuthal current, respectively, and the VEV 
of the radial component vanishes. These VEVs are periodic functions of the parameter a with 
the period equal to 1. So, if we present this parameter as (13. 7|) . with no being an integer 
number, then the VEVs are functions of qq alone. These functions are odd with respect to the 
reflection oq — ?• —oq- Both charge density and azimuthal current exhibit the jump structure at 
half-integer values for the ratio of the magnetic flux to the flux quantum. Simple expressions 
for the renormalized VEVs, Eqs. (I3.26P and (I3.28p . are obtained in the special case where 
the parameter q related to the planar angle deficit is an integer and the magnetic flux takes 
special values given by Eq. (]3.23p . In the general case of parameters a and q, we have derived 
two different representations for the renormalized VEVs. The first one is based on Eq. ()A.6P 
and the corresponding expressions for the charge density and azimuthal current have the forms 
()3.3ip and ()3.37p . The second representation is obtained by using the Abel-Plana summation 
formula and the corresponding expressions are given by Eqs. (I3.43p . (I3.46p . For a massless field 
the VEV of the charge density vanishes for points outside the magnetic vortex and the VEV 
of the azimuthal current behaves as r~^. For a massive field, for points near the vortex the 
VEVs behave as 1/r and for the charge density and the azimuthal current, respectively. 
At distances larger than the fermion Compton wavelength, the VEVs decay exponentially with 
the decay rate depending on the opening angle of the cone. The total charge induced by the 
magnetic vortex does not depend on the angle deficit of the conical space and is given by Eq. 
(I3.40p . In the special case of a magnetic flux in Minkowski spacetime, the formulae for the 
VEVs of the charge density and the azimuthal current reduce to Eqs. ()3.4ip . or equivalently, 
Eqs. (I3.47p . In Appendix |B] we show the equivalence of Eqs. (]3.47p to the expressions for the 
VEVs of charge density and azimuthal current known from the literature. 

The effects of a circular boundary on the VEV of the fermionic current are considered in 
Sect. [Hfor the exterior region. From the point of view of the physics in this region, the circular 
boundary can be considered as a simple model for the defect core. The mode sums of the charge 
density and the azimuthal current are given by Eqs. (14. 3 p and the radial component vanishes. In 
order to extract from the VEVs the contributions induced by the boundary, we have subtracted 
the boundary-free parts. Rotating the integration contours in the complex plane, we have derived 
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rapidly convergent integral representations for the boundary-induced contributions, Eqs. (j4.7p 
and (j4.12p . These formulae are simplified in the case of a massless field with expressions ()4.14p . 
In the exterior region, the total VEV of the fermionic current is a continuous function of the 
magnetic fiux. Note that unlike the boundary-free part, the boundary-induced part in the VEV 
of the charge density for a massless field is not zero. The parts in the VEVs induced by the 
boundary are periodic functions of the magnetic flux with the period equal to the flux quantum. 
These parts vanish for the special case of the magnetic flux corresponding to |ao| = 1/2. In 
the limit when the radius of the circle goes to zero and for |ao| < 1/2, for a massive fleld the 
boundary-induced contributions in the exterior region behave as a^'^", with qa defined by Eq. 
(j4.17p . For a massless field the corresponding asymptotics are given by Eqs. (|4.18p - (|4.20p . At 
large distances from the boundary and for a massive field, the contributions coming from the 
boundary decay exponentially [see Eqs. (I4.2ip ]. In the same limit and for a massless field, the 
boundary-induced VEV in the charge density decays as (a/r)'^'^". For the azimuthal current, 
the contribution induced by a circular boundary behaves as (a/r)^"^""*"^ when qa > 1/2, and like 
(a/r)^''" for qa < 1/2. Note that, when the circular boundary is present, the VEVs of physical 
quantities in the exterior region are uniquely defined by the boundary conditions and by the 
bulk geometry. This means that if we consider a non-trivial core model for both conical space 
and magnetic fiux with finite thickness b < a and with the line element (j2.ip in the region r > b, 
the results in the region outside the circular boundary will not be changed. 

The boundary-induced VEVs in the region inside a circular boundary are studied in Sect. 
El The corresponding mode sums for the charge density and the azimuthal current are given 
by Eq. (15. Sp . They contain series over the zeros of the function given by Eq. (15. 2p . For the 
summation of these series we have employed a variant of the generalized Abel-Plana formula. 
The latter allowed us to extract explicitly from the VEVs the parts corresponding to the conical 
space without boundaries and to present the contributions induced by the circle in terms of 
exponentially convergent integrals. In the interior region the boundary-induced parts in the 
renormalized VEVs of the charge density and the azimuthal current are given by Eqs. (15. lip 
and ()5.12p . For a massless fermionic field these formulae are reduced to Eqs. (j5.14p . For 
large values of the circle radius and for a massive field the boundary-induced contribution decay 
exponentially [Eqs. (|5.17p ]. In the case of a massless field, the VEVs decay as (r/a)^"?""^, for 
both charge density and azimuthal current. For points near the apex of the cone, the leading 
terms in the corresponding asymptotic expansions are given by Eqs. ()5.19p and ()5.20p . In 
particular, the boundary-induced parts vanish at the apex when |ao| < 1/2 — l/(2q) and diverge 
for |ao| > 1/2- l/(2g). 

The formulas for the VEV of the fermionic current are easily generalized for the case of a 
spinor field with quasiperiodic boundary condition (I2.28P along the azimuthal direction. This 
problem is reduced to the one we have considered by a gauge transformation. The corresponding 
expressions for the VEVs are obtained from those given above changing the definition of the 
parameter a by Eq. (j2.29p . 

The results obtained in the present paper can be applied for the evaluation of the VEV of the 
fermionic current in graphitic cones. Graphitic cones are obtained from graphene sheet if one 
or more sectors are excised. The opening angle of the cone is related to the number of sectors 
removed, N^, by the formula: cpo = 27r(l — Nc/Q) , with Nc = 1,2,... ,5 (for the electronic 
properties of graphitic cones see, e.g., [50] and references therein). All these angles have been 
observed in experiments [51] . The electronic band structure of graphene close to the Dirac points 
shows a conical dispersion -E'(k) = where k is the momentum measured relatively to the 

Dirac points and vp ^ 10^ cm/s represents the Fermi velocity which plays the role of a speed of 
light. Consequently, the long- wavelength description of the electronic states in graphene can be 
formulated in terms of the Dirac-like theory in (2+l)-dimensional spacetime. The corresponding 
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excitations are described by a pair of two-component spinors, corresponding to the two different 
triangular sublattices of tlie honeycomb lattice of graphene (see, for instance, |40]). In both cases 
of finite and truncated graphitic cones the corresponding 2-dimensional surface has a circular 
boundary. As the Dirac field lives on the cone surface, it is natural to impose bag boundary 
condition (|2.4p on the bounding circle which ensures the zero fermion flux through the edge of 
the cone. A more detailed investigation of the fermionic current in graphitic cones, based on 
the results of the present paper, will be presented elsewhere. 
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A Integral representations 

In this section we derive two integral representations for the function X{q, a^, z) defined by Eq. 
(I3.19p . In the first approach, we use the integral representation for the modified Bessel function 
1/3 .{z) (see formula 9.6.20 in Ref. [52J). To be allowed to replace the order of the integration 
and the summation over j, we integrate by parts the first term in this representation: 

^ , , sin(7r/3o) _„ z T , . sin(/3,,-?/) sin(7r/3,-) f°° , _,™h«-«« /a -.n 

Substituting (|A.ip into Eq. (|3.19p and interchanging the order of the summation and integration, 
we apply the formula |46j 

sin(/3jy) _ ^ ^y ir cos [(2/ + l)7r(ao - l/2g)] ^^^^ 
^ (3j qcos[7r{ao — l/2q)] ' 

for 2171 /q < y < (2/ + 2)71 /q. For the first term in the right-hand side of (jA.ip one has y = tt. 
For this term, when q = 21 with Z = 1, 2, . . ., the corresponding summation formula has the form 

^^^^Mil = {-1)1/'^- cos(g7rao) tan[7r(ao - l/2q)]. (A.3) 

Finally, for the series corresponding to the last term in Eq. (lA.ip we have 

^sm(vr/3,)e-'^^^ = J^^;"°'^\ (A.4) 
^ cosayqy) — cos(g7r) 

with the notation 

f{q,ao,y) = cos [gvr (1/2 - Qo)] cosh [(gao + g/2 - 1/2) y] 

- cos [qTT (1/2 + ao)] cosh [{qao - q/2 - 1/2) y] . (A. 5) 

Combining the formulae given above, we find the following integral representation for the 
series (13.19P : 



l{q,ao,z) = / dy- 

q TT Jq cosh(gy) - cos(Q'7r) 

+ - ^{-ly cos[27Tl{ao - l/2g)]e"'=°^(2"'/'?\ (A.6) 



1=1 
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with 2p < q < 2p + 2. In the case q = 2p, the term 



-(-l)''/'^sm(™), (A.7) 

should be added to the right-hand side of Eq. (|A.6p . Note that for I ^ q < 2, the last term on 
the right-hand side of Eq. (jA.Gh is absent. Formula ()A.6j) is simplified in the case q = 1: 

X(l,ao,.) = e- Tdye— (A.8) 
TT Jq cosh(y/2) 

In the special case (j3.23p with an integer q, the integral term in Eq. ()A.6P vanishes. For 

q = 2p + 1 one finds a, z) = {2/q) X][=o^^^^' e^™'^*-^'^'/'^^ , where, as before, the prime on the 
summation sign mens that the term with I = should be halved. For even values of q, by taking 
into account the additional term (|A.7p . we find I{q, a, z) = {2/q) Yl'i=o e^^°'^(^'^'/'') —e~^/q. Note 
that for an integer q, from definition (j3.19p one has T{q, a, z) = 2 Yl'^=o ^qn{z). Now, we can see 
that in the special case under consideration, formula ()A.6p coincides with Eq. (j3.24p . 

We can give an alternative integral representation of the series (j3.19p by using the Abel-Plana 
summation formula in the form (see \27\ H9] ) 

E /(- + 1/2) = r dz f{z) -^^dz ^^'"\;J}~''\ (A.9) 
n=o -^0 e +i 

First of all we write this formula in the form more appropriate for the application to Eq. (j3.19p . 
Let us consider the series 

oo 

Y^f(^\j + u\+ve,)= Y^Y.f{n+l/2 + 5\u + v\), (A.IO) 

j <5=±1 n=0 

where \u\ ^ 1/2, \v\ ^ 1/2. Applying formula ()A.9p . we get 

Y^f[\j + u\+vej)=J2 dzf{z + 5w)-tY, dz-^ Ji , ^ 

(A.ll) 

with the notation w = \u + v\. Introducing new integration variables, we present this formula 
in the form 

y2f{\j + n\+ve,) = 2 dzfiz)- dz [f (z) - f {- z)] 
, Jo Jo 



Now, deforming the integration contour, we write the integral along the half-line {i5w^ oo -|- i8w) 
in the complex plane z as the sum of the integrals along the segment {i6w,0) and along (0,oo). 
At this step we note that in the case 1/2 < u; < 1 the integrand has a pole at z = i5{w — 1/2). 
We exclude the poles by small semicircles in the right-half plane with radius tending to zero (see 
Fig. [6]). The sum of the integrals along the segments {iw^G) and {—iw,Q) cancels the second 
integral in the right-hand side of Eq. (IA.12P and we get the following result 

■ "'0 JO _5^_^^ ^ -r 
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where ^4 = for ^ |ii + f | ^ 1/2, and 

A = f{l/2-\u + v\)- f{\u + v\-l/2), 



(A.14) 



for 1/2 < |n + f I < 1. The term A comes from the contributions of the above mentioned poles 
to the integrals. 



iw 



-iw 



Figure 6: Integration contour in the complex plane z in the derivation of summation for- 
mula ()A.13p . 



We apply to the series in Eq. (13.190 formula (|A.13p with the function f[z) = Iqz{x) and 
with the parameters n = ao, v = —l/{2q). This leads to the following result 



2 

I{q, ao,x) = A{q, ao, x) + - dz Iz{x) 

Q Jo 

smh.{zTT)Kiz{x) 



-^Q Jo 



dz Re 



^2n{z+i\qao-l/2\)/g _|_ l 



with A{q,ao,x) = for \ao — l/2q\ ^ 1/2, and 



A{q,ao,x) = - sin[7r(|gao - 1/2| - q/2)]K\ga^^_i/2\-g/2{x), 



vr 



(A.15) 



(A.16) 



for 1/2 < |qo — l/2q| < 1. Note that in Eq. (jA.15p the function Kiz{x) is real and the real part 
in the integrand can be written explicitly by observing that 



Re(e 



U+IV 



+ 1) 



1 e " + cos V 



2 cosh u + cos V 
In the special case q = l, from (|A.15p we have 

j-OO 

Z{\,aQ,x) = A{l,aQ,x) + 2 I dz Iz{x) 

Jo 

smh.{zTr) Kiz{x) 



(A.17) 



4 f°° 
H — / dzRe 
71" Jo 



o27r(2-iao) _ \ 



(A.18) 



with A{l,ao,x) = — (2/7r) sin(7rao)-f^ao (x) in the case —1/2 < ao < and A{l,ao,x) = for 
^ ao ^ 1/2. Note that the last term in the right-hand side is an even function of ao. 
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B Comparison with the previous results in the Minkowski bulk 



The VEVs of the fermionic current induced by a magnetic flux in (2+l)-dimensional Minkowski 
spacetime have been previously investigated in Ref. |35J. The expressions for the VEVs of 
charge density and the azimuthal current derived in this paper have the form (in our notations) 

•0. ^^ _ _sin(7rao) r ^.JK^kr) 



{j {x))o,Ten = -em / dk 

71" Jm Vk^ - 



In this section, we show that these representations are equivalent to Eqs. (j3.47p . First we 
consider the charge density. As the first step, we use the integral representation for the square 
of the Macdonald function: 



\ r-e-^'"'-'K^,{z). (B.2) 



This formula is easily obtained from the integral representation of the product of two Macdonald 
functions given in Ref. [53j (page 439). Inserting Eq. (|B.2p into (jB.ip . after the change of the 
order of integrations, the integral over k is taken simply and for the charge density we obtain 
the result given in Eq. p.47p . 

In order to see the equivalence of the representations for the azimuthal current, we present 
the expression with the Macdonald functions in the form 

2 

Kl^{kr) - K^,_i{kr)K^,+i{kr) = j dxxKl^{kx). (B.3) 

After substituting this into Eq. (|B.ip . we use the integral representation (|B.2p . Then, we first 
take the integral over k and after that the integral over x. As a result, the integral representation 
(|3.47p for the azimuthal current is obtained. Hence, we have shown that, in the special case of 
Minkowski bulk, our results for the VEVs of the charge density and azimuthal current, given 
by Eq. (j3.47p . agree with those from the literature. Note that we have also derived alternative 
representations (13.4ip . For the numerical evaluation the latter are more convenient. 

C Contribution of the mode with j = —a 

When the parameter a is equal to a half-integer, the contribution of the mode with j = —a to 
the VEV of the fermionic current should be evaluated separately. Here we consider the region 
inside a circle with radius a. Similar to Eq. ()3.6p . the negative-energy eigenspinor for this mode 
has the form 

^'''-"^^"^A^ 1^ e^'?*/2cos(7r-7o) )' ^ ' 

where 70 is defined after Eq. (j3.6p . From boundary condition (|2.4p it follows that the eigenvalues 
of 7 are solutions of the equation 

m sin(7a) -|- 7 cos(7a) = 0. (C2) 

The positive roots of this equation we denote by 7^ = 7a, / = 1,2,.... From the normalization 
condition, for the coefficient in (jC.ip one has 

bl = [1 - sin(27a)/(27a)]-i . (C.3) 
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Using Eq. (jC.ip . for the contributions of the mode under consideration to the VEVs of the 
fermionic current we find: 



_ e 1 + /X [7; sin(27;r/a) - ^ cos(27;r/a)] /{aE) 

{J [x))j=_^ - 



ar<j)o fr[ 1 - sin(27,)/(27i) 



1-2, _ e 7; ^sin(27;r/a) + 7;cos(27;r/a) 

^''^^'=-" " ar^,^^{aEY 1 - sin(270/(270 ' ^ ^ ^ 

where {aE)"^ = '^f + and /i = ma. The part corresponding to the radial component vanishes. 
We assume the presence of a cutoff function. For the summation of the series in Eqs. (|C.4p . we 
use the Abel-Plana-type formula 



oo 



^ 1 - sm(27,)/(27l) l/l^+l " ' J„ fi^e=' + 1 

Eq. (jC.SP is obtained from more general summation formula given in [53] (see also j49] ) taking 
61 = and 62 = — For the functions /(z) corresponding to Eq. (|C.4p one has /(O) = 0. 
The contribution of the last term in Eq. (]C.5P to {j'^{x))j=-a is finite in the limit when the 
cutoff is removed. Noting that for both series in Eq. (IClill the function f{z) is an even function, 
we conclude that the last term in (jC.SP does not contribute to the both charge density and 
azimuthal current. For the remaining parts coming from the second term in the right-hand side 
of Eq. ()C.5p one has 

= — - / d7[l + M7)], 
Trrc^o Jo 



00 



0-^(x)),=_« = — d7[cos(27r) + /i(7)], (C.6) 

with /i(7) = mE^"^ [7sin(27r) — m cos(27r)] and i?^ = 7^ + m?'. These parts do not depend on 
the circle radius. Now, we can see that the integral with the function /i(7) is zero. The remained 
term in the charge density is subtracted by the renormalization and the remaining integral in 
the expression of azimuthal current is zero for r > 0. Hence, we conclude that the term with 
j = —a does not contribute to the renormalized VEV of the fermionic current. 
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